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Abstract. We summarize the main features of a class of anomalous (asymptotically flat, 
but non Schwarzschild-like) gravitational configurations in models of gravitating non-linear 
electrodynamics (G-NED) whose Lagrangian densities are defined as arbitrary functions of the 
two field invariants and constrained by several physical admissibility conditions. This class of 
models and their associated electrostatic spherically symmetric black hole (ESSBH) solutions 
are characterized by the behaviours of the Lagrangian densities around the vacuum and at the 
boundary of their domain of definition. 
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| 1. Preliminaries 

For several decades a great deal of attention has been focussed on G-NED as generalizations 
of the Reissner-Nordstrom (RN) solution of the Einstein-Maxwell field equations, whose more 
popular example is the gravitating Born-Infeld theory pQ. Some other G-NED models supporting 
asymptotically Schwarzschild-like solutions have been also considered in the literature [2]. In 
addition to such classes of models (and gravitational configurations) there are others leading 
to ESSBH solutions approaching asymptotic flatness slower than the Schwarzschild field In 
this work we shall analyze the full class of physically admissible G-NEDs supporting this kind 
of asymptotically anomalous solutions. 

The action determining the dynamics of a G-NED is defined as 



S = S G + S NE d = J d^xy 



R <P(X,Y) 



(1) 
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where ip(X,Y) is a (unspecified) function of the two field invariants X = —^F^F^ = 
E 2 — H 2 ,Y = —\F^ V F*^ U = 2E ■ H, which are built from the field strength tensor F^ v and its 

dual F* = \e^ va pF a ^ and where the electric E and magnetic H fields are defined in the usual 
way. The associated energy-momentum tensor in flat space is obtained as 

V = 2F m + ^F* v a ) - <p{X, Y)n, u . (2) 

Consistently with previous conventions we shall call physically admissible the models satisfying 
the minimal requirements of (i) definiteness, continuity and regularity of ip(X, Y) in their domain 



of definition of the X — Y plane, assumed to be open, connected and including the vacuum, (ii) 
parity invariance <p(X,Y) = <p(X,—Y), (iii) positive definite character of the energy, endorsed 
by the necessary and sufficient condition 



p > ( V^Tv 2 " + X) || + Y^ - <p(X, Y) > 0. (3) 

The associated flat-space field equations lead, for the electrostatic spherically symmetric (ESS) 
solutions (E(r) = E(r)r/r,H = 0), to the first-integral 



r 2^ 

dX 



__E{r) = Q, (4) 

where Q is an integration constant, identified as the electric charge. This equation allows the 
complete determination of the ESS solutions for admissible models. Assuming a field behaviour 
of the form E{r — > 0) ~ r p we obtain three possible cases: ultraviolet divergent (UVD) when 
p < — 1, for which the integral of the energy diverges around r = 0; Al when — 1 < p < 0, for 
which the field diverges at r = but the integral of energy converges there; A2 when p = 
(and E(r —5-0) ~ a — br a with a > 0) for which the central field is finite and the integral 
of energy converges there. Concerning the asymptotic behaviour of the solutions, assumed to 
be of the form E(r —> oo) ~ r q , we are lead to four possible cases: infrared divergent (IRD) 
when — 1 < q < 0, for which the fields vanish as r — > oo but the integral of energy diverges 
there; Bl (—2 < q < —1), B2 (q = —2) and B3 (q < —2) cases, for which the integral 
of energy converges asymptotically. The twelve combinations of these central and asymptotic 
behaviours exhaust the full class of admissible NEDs. Excluding the models leading to UVD 
and IRD behaviours we obtain six classes exhausting the set of admissible models supporting 
finite-energy ESS solutions. The minimal coupling to gravity of the nine classes of models 
obtained by excluding the IRD behaviour was extensively analyzed in Ref.jlj. This analysis 
determines and classifies the possible G-ESS configurations, which are reduced to extreme BHs 
(EBH), one or two-horizons BHs, extreme and non-extreme black points and naked singularities 
(NS), depending on the values of Q and the ADM mass coming from the integration of the 
field equations. All these configurations behave asymptotically as the Schwarzschild field. As 
we shall see, the three classes of models supporting ESS solutions with IRD behaviour lead to 
asymptotically anomalous gravitational configurations with similar structures at finite r. 

2. Asymptotically anomalous gravitational configurations 

Let us write the general expression of the G-ESS solutions associated to the action (pQ). The 
symmetry of the source Tq = leads to the static spherically symmetric line element [1] 

ds 2 = \{r)dt 2 - \- l {r)dr 2 - r 2 dtt 2 , (5) 

where dQ 2 = d9 2 + sin 2 0d(/? 2 . With the line element © both the first-integral and 
X = E 2 (r,Q) have the same expression in flat space and in the gravitational context, allowing 
an immediate translation of all results obtained in the former case to the latter one. 

Let us analyze three classes of admissible models supporting asymptotically IRD solutions 
with central field behaviours AI, A2 or UVD. Now the integral of energy diverges at r — > oo as 
r q+1 , (— 1 < q < 0) but converges at the center in the AI and A2 cases, allowing the definition 
of the interior integral of energy Si n (r, Q) = J*q R 2 Tq(R, Q)dR which, owing to the admissibility 
conditions, is a monotonically increasing and convex function of r. The general solution 
of the Einstein equations for the metric reads in these cases 

X(r,Q,C) = 1 + — , (6) 



Figure 1. : Interior integral of energy for the different families of NED models with 
asymptotically IRD field behaviours, cut by the beam of straight lines (r + C)/2. 



where C is an integration constant. The radii of the horizons (zeroes of ([6])) are obtained as 
solutions of = Ei n (r, Q) and thus are given by the intersections between the curves £» n (r, Q) 
and the beam of straight lines (r + C)/2 corresponding to different values of C and Q (see Fig.l). 
Through this procedure we are led to the following structures for the configurations: 



2.1. Al-IRD 

In this (r, Q) has a divergent slope as r ~ and the above method leads to five classes 

of configurations. In particular there is a value C = C ex t r (Q), corresponding to the straight line 
tangent to the curve £j n (r, Q), defining the radius r^extriQ) of a EBH configuration (for which 
the horizon is degenerate, i.e. A(r) = \'{r)\ r=rh = 0). These values are obtained as 

8nrl extr (Q)T$(r hextr ,Q) = I ; C extr {Q) = 2e in (r hextr (Q),Q) -r hextr (Q). (7) 

As easily seen in Fig.l, for C > C ex t r (Q) there are not cut points and the associated 
configurations are NS, while for < C < C ex t r {Q) they are two- horizons BHs. On the other 
hand, for C < there is always a single cut point corresponding to BHs with a single (non- 
extreme) horizon. In the critical case (C = 0) the metric diverges to —00 around r = as 

A(r) = 1 - , 32 ^ 9 N rP + • • ■, (8) 

and increases monotonically to A(r — > 00) = 1 exhibiting a single- horizon as well. Consequently 
the central singularity is timelike for C > and spacelike for C < 0. As C — > in Eq.© 
the metric function converges to the critical one for r > but differs at the center by 5-like 
terms. At large r, for all these IRD cases the asymptotic behaviour of the metric ^ approaches 
flatness as A(r — > 00) — 1 ~ —r q (for — 1 < q < 0, and as ~ ln(r)/r for q = — 1) and is not 
Schwarzschild-like. Thus the ADM mass is not defined in the asymptotically "anomalous" cases. 



2.2. A 2- IRD 

From Fig.l we see that similar configurations as in case Al are obtained for while for the 
critical configuration, C = 0, the metric at the center, which behaves as 

32vr6Q 



A(r 0, C = 0, Q, a, b, a) ~ 1 - 16vrQa + 



-r + x r 



(9) 



(x- integration constant) is finite. Its sign there depends on the quantity 16irQa. In the A2a 
case the critical metric takes the value at the center A(0) = 1 — 16nQa > and its slope there is 



dr 



32irbQa 



r^o (a+l)(2-ff 



-r^ 1 + 2xr. 



(10) 



This implies that the slope of the metric at the center for the critical configuration is non-negative 
and depends on a, becoming divergent for a < 1, vanishing for a > 1 and taking a finite value 
for a = 1. As a consequence the critical metric, which is always monotonically increasing, leads 
to timelike NS. On the other hand, in case A2b corresponding to A(0) = 1 — 16irQa < the 
critical configuration exhibits a single-horizon BH with a spacelike singularity at the center. The 
case A2c with A(0) = 1 — lQirQa = leads to extreme black points [5]. 



2.3. UVD-IRD 

In this (r, Q) cannot be defined, but now the Einstein equations can be integrated as 

where e(r,Q,T) = Air j r 2 T$(r,Q) + T and the arbitrary constant T has been absorbed in 
the constant D. This function is monotonically increasing and convex, exhibiting a vertical 
asymptote at r = and diverging with vanishing slope as r — > 00. Following the same procedure 
as in the previous cases, the radii of the horizons are obtained from the cutting points of the 
function e(r, Q,0) with the beam of straight lines (r + D)/2, corresponding to different values 
of D. There is again a tangent line to the curve e(r, Q,0), determining the radius of an EBH 
configuration through the same equation ([7]). For D > Dhextriff) the configurations are NS, 
while two-horizons BHs appears for D < Di iextr {Q), as can be seen from Fig.l. 



3. Discussion 

The analysis of the asymptotically anomalous class of NED models considered here, together 
with the one of the asymptotically Schwarzschild-like class studied in exhausts the family 
of G-ESS solutions of physically admissible NEDs. For these families, aside from NS, EBHs 
and two-horizon BHs already present in the RN case, there are in addition single-horizon BHs 
and (in case A2) extreme and non-extreme black points [5]. Thermodynamics of the different 
gravitating ESS structures outlined here, making use of the fact that NED models satisfy both 
the zeroth and first laws of BH thermodynamics [6j, is currently in progress. 

Let us finally point out that the present analysis can be extended to non-abelian gauge fields. 
In this case the lagrangian density is still given by Eq.([T]), where the field invariants are obtained 
by using the usual prescription in the calculation of the traces over the gauge group generators, 
leading to X = — ^F^ ua F^ u = E 2 a -H 2 a ,Y = -±F^ ua F^ u = 2E a -H a , a = N (however, other 

prescriptions are possible). This problem can be reduced, for ESS fields {E a (r) ^ 0, H a = 0), to 
the abelian one in a correspondence between ESS fields of both theories associated to the same 
form of the Lagrangian density, solving the gravitating problem for the former in terms of the 
solution of the latter. This extension will appear in a forthcoming publication. 
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